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Abstract 

The first measurement on the antiproton to proton ratio by the AMS-02 col¬ 
laboration agrees with the expectation from conventional cosmic-ray secondaries in 
the kinetic energy range ~ 20 — 100 GeV, which can be turned into stringent upper 
limits on the annihilation cross sections of dark matter (DM) above ~ 300 GeV. Us¬ 
ing the GALPROP code, we derive the upper limits in various propagation models 
and DM profiles. We show that in the “conventional” propagation model, for the 
qq, bb, and W + W~ final states, the constraints could be more stringent than that 
derived from the recent Ferm-LAT gamma-ray data on the dwarf spheroidal satellite 
galaxies. Making use of the typical minimal, median and maximal models obtained 
from a previous GALPROP based global fit to the preliminary AMS-02 data, we 
show that the variation of the upper limits is about a factor of five. The possibility 
of DM contributions to the low and very high energy p/p data is discussed. 
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Dark matter (DM) is known to contribute to 26.8% of the total energy density of the 
Universe [I]. However, its particle nature remains largely unknown. If the DM particles 
can annihilate or decay into the standard model (SM) final states, they will contribute 
to new primary sources of cosmic-ray particles, and result in significantly changes in the 
spectra of the cosmic-ray antiparticles such as positrons and antiprotons, as these species 
are assumed to be secondaries in the conventional cosmic-ray propagation theory. An¬ 
tiprotons are highly expected from DM annihilation in many DM models, and unlikely to 
be generated from the nearby pulsars. Compared with cosmic-ray electrons/positrons, the 
cosmic-ray antiprotons loss much less energy during propagation, and can travel through 
longer distance in the Galaxy, which makes the antiproton flux more sensitive to the 
uncertainties in the propagation parameters and the DM profiles. 

Recently, the AMS-02 collaboration has released the first preliminary result of the 
cosmic-ray antiproton to proton flux ratio ft/p [2j. The measured kinetic energies of the 
antiprotons have been extended to ~ 450 GeV. Although the spectrum of ft/p at high 
energies above 100 GeV tend to be relatively flat, within uncertainties the AMS-02 data 
are consistent with the background of secondary antiprotons, which can be used to set 
stringent upper limits on the dark matter (DM) annihilation cross sections, especially for 
high mass DM particles. The constraints on the DM properties from antiprotons have been 
investigated previously before AMS-02 ( see eg. dam la ©d )• In this work, we explore the 
significance of the new AMS-02 ft/p data on constraining the annihilation cross sections 
of the DM particles in various propagation models and DM profiles. Four representative 
background models are considered with four different DM profiles. We derive the upper 
limits using the GALPROP code and show that in the “conventional ” propagation model 
with Einasto DM profile, the constraints could be more stringent than that derived from 
the Ferm-LAT gamma-ray data on the dwarf spheroidal satellite galaxies. Making use of 
the typical minimal, median and maximal models obtained from a previous GALPROP 
based global fit to the preliminary AMS-02 data, we show that the uncertainties on the 
upper limits is around a factor of five. 

We start with a briefly overview on the main features of the cosmic-ray propaga¬ 
tion within the Galaxy. The Galactic halo within which the diffusion processes occur is 
parametrized by a cylinder with radius Rh = 20 kpc and half-height — 1 — 20 kpc. The 
diffusion equation for the cosmic-ray charged particles reads (see e.g. 0) 
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where i()(r,p,t) is the number density per unit of total particle momentum. For steady- 
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state diffusion, it is assumed that dijj/dt = 0. The number densities of cosmic-ray particles 
are vanishing at the boundary of the halo, i.e., ip(Rh, z,p ) = ±Z h ,p) = 0. The energy 
dependent spatial diffusion coefficient D xx is parametrized as D xx = /3D 0 ( p/poY , where 
p = pj(Ze) is the rigidity of the cosmic-ray particle with electric charge Ze. The power 
spectral index S can have different values 5 = Sim when p is below (above) a reference 
rigidity p 0 - The coefficient D 0 is a normalization constant, and f3 = v/c is the velocity of 
the cosmic-ray particle. The convection term in the diffusion equation is related to the 
drift of cosmic-ray particles from the Galactic disc due to the Galactic wind. The diffusion 
in momentum space is described by the reacceleration parameter D pp which is related to 
the velocity of disturbances in the hydro dynamical plasma, and described by the Alfven 
speed V a [8]. In Eq. (JT|) , the momentum loss rate is denoted by p which could be due 
to ionization in the interstellar medium neutral matter, Coulomb scattering off thermal 
electrons in ionized plasma, bremsstrahlung, synchrotron radiation, and inverse Compton 
scattering, etc.. The parameter Tf(r r ) is the time scale for fragmentation (radioactive 
decay) of the cosmic-ray nuclei as they interact with interstellar hydrogen and helium. 

The spectrum of a primary source term for a cosmic-ray nucleus A is assumed to have 
a broken power low behaviour dqA(p)/dp oc (p/p^ s ) 7A with 7^4 = 7 ^ 1 ( 7 , 42 ) for the nucleus 
rigidity p below (above) a reference rigidity Pa s ■ The spatial distribution of the primary 
sources is taken from Ref. [9|. Secondary antiprotons are created dominantly from inelastic 
pp- and pA-collisions with the interstellar gas. The corresponding source term reads 

,Vi(p,p') 
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where n, is the number density of the interstellar hydrogen (helium), n p is the number 
density of primary cosmic-ray proton per total momentum, and d<Ji(p,p')/dp' is the dif¬ 
ferential cross section for p + H(He) —> p + X. In calculating the antiprotons, inelastic 
scattering to produce “tertiary” antiprotons should be taken into account. 

The primary source term from the annihilation of Majorana DM particles has the 
following form 
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where (crv) is the velocity-averaged DM annihilation cross section multiplied by DM rel¬ 
ative velocity (referred to as cross section). p(r) is the DM energy density distribution 
function, and dN^ x ’/dp is the injection energy spectrum of antiprotons from DM an¬ 
nihilating into SM final states through all possible intermediate states X with r)x the 
corresponding branching fractions. 

The interstellar flux of the cosmic-ray particle is related to its density function as 
= v%l}(r,p)/(An). At the top of the atmosphere (TOA) of the Earth, the fluxes of 
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cosmic-rays are affected by solar winds and the helioshperic magnetic field. This effect is 
taken into account using the force-field approximation which involves a parameter </>, the 
so called Fisk potential [TO] . In this work, we shall take </> = 550 MV in the numerical 
analysis. 

We shall solve the diffusion equation of Eq. ([T]) using the publicly available code GAL- 
PROP v54 [11, [12] 13] QH [IS] which utilizes realistic astronomical information on the 
distribution of interstellar gas and other data as input, and considers various kinds of 
observables in a self-consistent way. Other approaches based on simplified assumptions on 
the Galactic gas distribution which allow for fast analytic solutions can be found in Refs. 

nsn. 

We first consider the so-called “conventional” diffusive re-acceleration (DR) model 
PIII5] which is commonly adopted by the current experimental collaborations such as 
PAMELA [21, 22| [23] and Fermi-LAT [24] [25] as a benchmark model for the astrophysical 
backgrounds. It is useful to consider this model as a reference model to understand how 
the DM properties could be constrained by the AMS-02 data. Then we consider three 
representative propagation models selected from a large sample of models obtained from 
a global Bayesian MCMC fit to the preliminary AMS-02 proton and B/C data using the 
GALPROP code [26]. They are selected to represent the typically minimal (MIN), median 
(MED) and maximal (MAX) antiproton fluxes within 95% CL, corresponding to the region 
enveloping 95% of the MCMC samples with highest likelihoods in a six-dimensional pa¬ 
rameter space. For simple Gaussian posterior distributions, it corresponds to Ay 2 = 12.59 
for six degrees of freedom. The parameters in the four models are summarized in Tab. Q] 

Note that the “MIN”, “MED” and “MAX” models used in this work are different 
from and complementary to that given in Ref. [27] in several ways: i) they are obtained 
using the fully numerical GALRPOP code which is consistent with the analysis framework 
of this work while that in Ref. 123 are based on the two-zone diffusion model with 
the assumption of simplified galactic geometry and uniform interstellar gas distributions, 
ii) they correspond to the DR propagation model while that in Ref. [27] the models 
with diffusive-reacceleration plus a constant convection velocity V c are considered. In 
the GALPROP approach, it is known that the spectral shape of the B/C flux ratio is 
better reproduced without including the convection term; iii) the models in this work are 
based on the global fit to the preliminary AMS-02 data of proton flux and B/C flux ratio, 
while the models in Ref. [27] are based on the much older HEAO-3 data on B/C ratio 
with significantly lower precision. Consequently, the uncertainties in the prediction for 
the DM induced antiproton fluxes are much larger; iv) the Bayesian analysis was used in 
deriving these models, while in Ref. m the authors performed a frequentist analysis. The 
credible intervals obtained from the Bayesian analysis can be different from that using the 


4 


confidence levels based on a given value of Ay 2 . 

The predicted proton flux and the B/C flux ratio in these models are shown in Fig. [T] 
together with the latest AMS-02 data [2] • The figures show an overall agreement with the 
current data in these models. In the “conventional” model, the predicted B/C ratio is a 
little higher for the kinetic energy below ~ 10 GeV/n, but are consistent with the B/C 
data in the higher energies. Note that the theoretical predictions from Ref. [26] is based 
on the analysis of the previous preliminary AMS-02 results announced at the conference 
ICRC(2013) [28], [29] . So far the AMS-02 data on the B/C flux ratio is still unpublished. 
We will update the whole analysis in Ref. [26] after they get published, in a future work. 
The predictions for the background of the p/p flux ratio in these models are shown in 
Fig. [2] The “MIN”, “MED” and “MAX” models are highly degenerate in the background 
p/p ratio. Compared with these models, the “conventional” model predicts more low 
energy antiprotons but at high energies above ~ 500 GeV, the predicted antiprotons are 
much less. In all the four DR propagation models, below ~ 10 GeV the GALPROP 
based calculations underpredict the p/p flux ratio by ~ 40%, which is a known issue. 
The agreement with the low energy p data can be improved by introducing breaks in 
diffusion coefficients ffl. “fresh” nuclei component EU or a DM contribution [3]. Note 
that at energies below ~ 10 GeV, the uncertainties due to the solar modulation and the 
propagation parameters increase significantly. The predictions for low energy p/p ratio can 
be easily modified by introducing an independent Fisk potential 0 for p and an energy- 
dependent overall normalization factor as discussed in Ref. [32]. In this work, we instead 
use these DR models to derived very conservative upper limits on the annihilation cross 
sections of light DM particles. Note however that in the DR propagation models, the 
background predictions agree with the AMS-02 data well at higher energies ~ 10 — 100 
GeV, which can be turned into stringent constraints on the nature of heavy DM particles. 


model 

A(kpc) 

Z h {kpc) 

Do 

Po 

Si/5 2 

K(km/s) 

Ps 

7pi/ 7 P 2 

Conventional 

20 

4.0 

5.75 

4.0 

0.34/0.34 

36.0 

9.0 

1.82/2.36 

MIN 

20 

1.8 

3.53 

4.0 

0.3/0.3 

42.7 

10.0 

1.75/2.44 

MED 

20 

3.2 

6.50 

4.0 

0.29/0.29 

44.8 

10.0 

1.79/2.45 

MAX 

20 

6.0 

10.6 

4.0 

0.29/0.29 

43.4 

10.0 

1.81/2.46 


TAB. 1: Parameters in the propagation models “Conventional” pans], “MIN”, “MED” 
and “MAX” models from Ref. [26]. Dq is in units of 10 28 cm 2 • s -1 , the break rigidities po 
and p s are in units of GV. 

The flux cosmic-ray antiprotons from DM annihilation depend also significantly on the 
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FIG. 1: Predictions for the proton flux (left) and the B/C flux ratio (right) in the four 
propagation models listed in Tab. [H The latest data of proton flux from AMS-02 [2] and 
PAMELA [33, US] are shown. 


choice of DM halo profile. N-body simulations suggest a universal form of the DM profile 


P(r) 


( 7 ( 1 + Po/ps) 0 ^ (/3 7)/l 

p& V r ©y \ l + ( r lr Q ) a ) 


(4) 


where p Q ~ 0.43 GeV cm -3 is the local DM energy density [35j. The values of the pa¬ 
rameters a, (3, 7 and r s for the Navarfro-Frenk-White (NFW) profile [36], the isothermal 
profile [37 and the Moore profile [ 351 [33 j are summarized in Tab. El An other widely 



a 

P 

7 

A (kpc) 

NFW 

1.0 

3.0 

1.0 

20 

Isothermal 

2.0 

2.0 

0 

3.5 

Moore 

1.5 

3.0 

1.5 

28.0 


TAB. 2: Values of parameters a, f3, 7 and r s for three DM halo models, NFW [36] , 
Isothermal [37] , and Moore (381 09] • 


adopted DM profile is the Einasto profile m 


p(r) = p 0 exp 



rp&E _ r a E 

' A o 


• S 


1 


(5) 


with cxe ~ 0.17 and r s « 20 kpc. 

We consider three reference DM annihilation channels XX ~t XX where = qq , bb 
and W + W~. The energy spectra of these channels are similar at high energies. The main 
difference is in the average number of total antiprotons N x per DM annihilation of each 
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FIG. 2: Predictions for the p/p ratio from the four propagation models list in Tab. [0 
The data from AMS-02 [2] and PAMELA [33] are shown. 


channel. For a DM particle mass m x = 500 GeV, the values of N x for typical final states 
are N qq = 2.97 (q = u,d ), N bb = 2.66, and Nww = 1-42. The injection spectra dN^/dp 
from DM annihilation are calculated using the numerical package PYTHIA v8.175 0B- 
in which the long-lived particles such as neutron and are allowed to decay and the 
final state interaction are taken into account. Since PYTHIA v8.15 the polarization and 
correlation of final states in r-decays has been taken into account [42]. 

In this work, we shall first derive the upper limits on the DM annihilation cross section 
as a function of DM particle mass, using the frequentist y 2 -analyses. The expression of 
y 2 is defined as 


(f th fexp\2 

2 _ Ut h ) 

X / . 9 5 


( 6 ) 


where f/ h are the theoretical predictions. f/ xp and cr* are the central values and errors of 
experimental data, respectively. The index i runs over all the available data points. For 
a given DM particle mass, we first calculate the minimal value Xmm °f the y 2 -function, 
and then derive the one-side 95% CL upper limits on the annihilation cross section, corre¬ 
sponding to Ay 2 = 3.84 for one parameter. All of the 30 data points of the AMS-02 p/p 
data are included in calculating the limits. 

In Fig. [3] we show the obtained upper limits on the cross sections for DM particle 
annihilation into bb final states from the AMS-02 p/p data in the “conventional”, “MED”, 
“MIN” and “MAX” propagation models. Four different DM profiles NFW [36], Isother¬ 
mal m, Einasto [40] and Moore [381 39] are considered. As can be seen, the upper limits 
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FIG. 3: Upper limits on the cross sections for DM particle annihilation into bb final 
states from the AMS-02 p/p data in the “conventional” (upper left), “MED” (upper 
right), “MIN” (lower left) and “MAX” (lower right) propagation models. Four DM pro¬ 
files NFW [36], Isothermal [37j , Einasto jlO] and Moore [38, [39] are considered. The 
upper limits from the Fermi-LAT 6-year gamma-ray data of the dwarf spheroidal satellite 
galaxies of the Milky Way are also shown [43]. The horizontal line indicates the typical 
thermal annihilation cross section (crv) = 3 x 10 -26 cm 3 s -1 . 


as a function of m x show some smooth structure for all the final states and DM profiles. 
The limits tend to be relatively stronger at m x & 300 GeV, which is related to the fact 
that the background predictions agree with the data well at the antiproton energy range 
~ 20 — 100 GeV. For a comparison, the upper limits from the Fermi-LAT 6-year gamma- 
ray data of the dwarf spheroidal satellite galaxies of the Milky Way are also shown [43] . 
In the “conventional” model, the upper limits from the AMS-02 p/p data are found to be 
compatible with that derived from the Fermi-LAT gamma-ray data for m x > 300 GeV. 
This observation holds for most of the DM profiles. In the “MED” model, the constraints 
are relatively weaker, which is related to the under prediction of low energy antiprotons 
in this model and the limits are more conservative. For an estimation of the uncertainties 
due to the propagation models, from the “MIN” model to the “MAX” model, we find that 

































FIG. 4: The same as Fig. [3j but for DM annihilation into W + W final states. 


the variation of the upper limits is within about a factor of five. 

For the W + W~ final states, the results are shown in Fig. [4] The constraints from 
AMS-02 p/p data turn out to be more stringent than that from the Fermi-LAT gamma- 
ray data for all the four DM profiles in the “convention” model when the DM particle mass 
is above ~ 300 GeV. Again we find that the variation of the upper limits from the “MIN” 
to the “MAX” model is within a factor of five. The result for the qq final states is shown in 
Fig. [5j Similar to the case of W + W~ final states, the constraints from p/p data are more 
stringent at about ~ 300 GeV. Compared with the case of W + W~ and bb, the constraints 
on the qq final states are the most stringent. For all the three final states, we find that 
the allowed DM annihilation cross section is below the typical thermal cross section for 
m x < 300 GeV in the conventional propagation model with Einasto profile, which shows 
that the AMS-02 p/p data can impose stringent constraints on DM candidates of weakly 
interacting massive particles. 

As can be seen from Fig. [21 compared with the AMS-02 data the GALPROP DR 
models predict fewer antiprotons at low (< 10 GeV) and very high (> 100 GeV ) energies. 
Without a robust estimation of the theoretical uncertainties, it is premature to claim any 
excesses in the p/p data. Nevertheless, we consider what would be the implications for 
DM if such a trend in the observations is confirmed by future analyses. The low energy 
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FIG. 5: The same as Fig. [3l but for DM annihilation into qq final states. 


data would allow for a non-vanishing DM annihilation cross section. For instance, in the 
“conventional” propagation model, for rn x =10.1, 35.0 and 75.8 GeV, the best-fit values 
are (av) = 3.6 x 10~ 27 , 1.14 x 10~ 26 , and 2.79 x 10~ 26 cm 3 s -1 , respectively, if the DM 
profile is Einasto, and the DM particles annihilate dominantly into bb final states. If both 
m x and (av) are allowed to vary freely, the best-fit DM particle masses and annihilation 
cross sections are 

m x = 58.5 (35.0) GeV and (av) = 2.16 (0.86) x 10" 26 cmV 1 

for DM annihilating into bb ( qq ) final states. In Fig. [6l we show the calculated spectra of 
p/p flux ratio from the best-fit DM particle masses and cross sections. The figure shows 
that the low energy p/p data are well reproduced by including such a DM contribution, 
except for the data point with kinetic energy below 1 GeV. 

As shown in Fig. [21 the spectrum of the AMS-02 p/p ratio tends to be flat toward 
high energies above ~ 100 GeV. This trend, if confirmed by the future AMS-02 data, 
is not expected from the secondary production of antiprotons, and raises the interesting 
question whether this would leave some room for a heavy DM contribution, similar to 
the case of the AMS-02 positron fraction jUJ EH EO :|8]. To explore this possibility, 
we perform an other fit using the p/p ratio data above 20 GeV (15 data points in total) 
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FIG. 6: Left) Spectrum of p/p flux ratio from DM annihilating into bb final states with 
m x = 58.5 GeV and (av) = 2.16 x 1CD 26 cm 3 s _1 obtained from a fit to the whole AMS- 
02 p/p data [2]. The “conventional” background model and the Einasto DM profile are 
assumed. Right) The same as left, but for the fit with qq final state with the best-fit values 
m x = 35 GeV and (av) = 0.86 x 10 -26 cm 3 s _1 . 


in order to avoid the theoretical uncertainties in the low energy region. The obtained 
Xmin as a function of m x for the bb, qq and W + W~ final states in the “conventional” 
propagation model with Einasto DM profile are shown in Fig. [TJ One can see that for 
the three final states the values of Xmm decrease almost monotonically from ~ 21 to 
~ 5 with an increasing DM particles mass from 100 GeV to 10 TeV, but the y 2 -curves 
become gradually flat toward high DM masses. Only for the W + W~ channel, there exists 
a shallow local minimal at around 6.5 TeV with low statistical significance. From the 
X 2 -curves, one can see that the DM particles mass is restricted to be above ~ 2 TeV at 
2cr. For an illustration purpose, we show in Fig. [7] the predictions for the p/p ratio in 
the “conventional” background model with a DM contribution. The DM particles masses 
and annihilation cross sections chosen to be m x = 6.5 TeV, (av) = 1.9 x lCV 24 cm 3 s -1 for 
W + W~, m x = 10.9 TeV, (av) = 3.4 x 10~ 24 cm 3 s _1 for bb channel, and m x = 10.9 TeV 
and (av) = 3.3 x 10~ 24 cm 3 s _1 for qq channel. Note that these values are not from the best- 
fit values. We conclude that introducing a DM contribution can improve the agreement 
with the AMS-02 p/p data with kinetic energy above 100 GeV, but the statistics is not 
high enough to determine the DM properties such as its mass and interaction strength. 
If the DM particle mass is indeed at 0(10) TeV scale, next generation precision cosmic- 
ray detection experiments are needed. As can be see in Fig. O the possible “excess” is 
located at the kinetic energy range 100 — 450 GeV where the secondary backgrounds from 
the four propagation models are similar. However, beyond ~ 450 GeV, the p/p from the 
“conventional” model drops quicker than that in the other propagation models. The future 
high energy antiproton data will be very important not only in probing DM but also in 
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constraining the background models. 

In summary, we have explored the significance of the first AMS-02 p/p data on con¬ 
straining the annihilation cross sections of the DM particles in various propagation models 
and DM profiles. Four representative background models have been considered with four 
different DM profiles. We have derived the upper limits using the GALPROP code and 
shown that in the “conventional ” propagation model with Einasto DM profile, the con¬ 
straints can be more stringent than that derived from the Ferm-LAT gamma-ray data on 
the dwarf spheroidal satellite galaxies. Making use of the typical minimal, median and 
maximal models obtained from a previous global fit, we have shown that the uncertainties 
on the upper limits is around a factor of five. The future more precise AMS-02 data can 
help to reduce the uncertainties in the derived upper limits. 

Note added: As we were completing this study, Ref. [32] appeared on the arXiv, which 
addresses some of the same problems as discussed here. Although the conclusions are 
similar, the analysis in this work is based on the fully numerical GALPROP code, while 
that in [32] is based on the two-zone diffusion model with (semi)-analytical approach. 
The two methods are quite complementary to each other. Compared with Ref [32], the 
upper limits obtained in this work are weaker at DM particle mass below ~ 100 GeV, but 
stronger for heavy DM particles above ~ 500 GeV for typical DM particles annihilating 
into bb final states. Similar discussions on the DM matter contributions can be found in 
Refs. [49] [50]. 
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FIG. 7: (Upper left) values of Xmin as a function of DM particle mass m x from a fit 
to the AMS-02 p/p data ( with kinetic energy above 20 GeV ) in the “conventional” 
propagation model pans] with the DM profile fixed to Einasto ffl- Three annihilation 
channels bb, qq and W + W~ are considered. (Upper right) predicted p/p ratio in the case of 
background (“conventional” model) plus a DM contribution with m x = 6.5 TeV, (av) = 
1.9 x 10 _24 cm 3 s“ 1 , and annihilation final states W + W~. The flux ratio of antiproton 
from DM to the proton from the background Pdm/pbg is shown as the dashed line. The 
data from AMS-02 [2] and PAMELA {33j are also shown. (Lower left) the same as the 
upper right, but for the bb channel with m x = 10.9 TeV and (av) = 3.4 x 10” 24 cm 3 s _1 . 
(Lower right) the same as the upper right, but for the qq channel with m x = 10.9 TeV 
and (av) = 3.3 x 10~ 24 cm 3 s _1 . 
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